Dzyaloshinskii-Moriya interactions in valence bond systems 
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We investigate the effect of Dzyaloshinskii-Moriya interactions on the low temperature magnetic suscepti- 
bility for a system whose low energy physics is dominated by short-range valence bonds (singlets). Our general 
perturbative approach is applied to specific models expected to be in this class, including the Shastry-Sutherland 
model of the spin-dimer compound SrCu2(B03)2 and the antiferromagnetic Heisenberg model of the recently 
discovered 5=1/2 kagome compound ZnCua(OH)eCl2. The central result is that a short-ranged valence bond 
phase, when perturbed with Dzyaloshinskii-Moriya interactions, will remain time-reversal symmetric in the ab- 
sence of a magnetic field but the susceptibility will be nonzero in the T — > limit. Applied to ZnCu3(OH)eCl2, 
this model provides an avenue for reconciling experimental results, such as the lack of magnetic order and lack 
of any sign of a spin gap, with known theoretical facts about the kagome Heisenberg antiferromagnet. 



I. INTRODUCTION 

A question at the heart of frustrated magnetism is what 
happens to an antiferromagnet when conventional magnetic 
states, such as the Neel state, are destabilized by frustrating in- 
teractions or geometry. Classically, a system commonly finds 
itself in a degenerate manifold of configurations which share 
a local constraint, and whose size grows exponentially with 
the system 1 . For quantum spins, this would violate the third 
law of thermodynamics. Instead, it has been suggested that 
the system will enter one of a number of exotic non-magnetic 
phases 2 ' 3 ' 4 , with orders that do not have classical analogs. 
While some of these exotic quantum phases have now been 
seen in toy model s 5 ' 6 ' 7 ' 8 ' 9,10 , settling the question of their exis- 
tence as a matter of principle, an unambiguous experimental 
observation of such phases in actual materials is still lacking. 
Therefore, the issue of what experimental signatures are in- 
dicative or at least suggestive of such phases is a question of 
considerable interest. The most basic feature is the absence 
of magnetic order even when the temperature is much lower 
than the antiferromagnetic coupling J. 

Within a non-magnetic phase, one may visualize the system 
wave function as a superposition of valence bond states, i. e. 
states where every spin is paired with another spin to form a 
singlet or valence bond. The valence bond states form a highly 
overcomplete basis for the 5 — subspace 11 so the represen- 
tation need not be unique. However, it is convenient to distin- 
guish between phases in which the wave function can be ex- 
pressed mainly in terms of short-range valence bond states 12 , 
where the singlet pairing is always between nearby spins, and 
phases where the description always involves important con- 
tributions from valence bonds of all lengths. In the latter case, 
the most famous example of which is the original formulation 
of the resonating valence bond (RVB) liquid 213 , the spin-spin 
correlation decays algebraically and the magnetic excitations 
are gapless. In the former case, which includes the valence 
bond solids 3 and the short-range Z2 and U(l) RVB liquids^! 4 ., 
spin-spin correlations decay exponentially. While, in princi- 
ple, the short-ranged nature of equal time correlations does not 
necessarily imply a spectral gap 15 (and hence it might be pos- 
sible to conceive of a phase with exponential spin-spin cor- 
relations and gapless magnetic excitations), the most familiar 
examples of short-range valence bond phases, including the 



three just mentioned, have a spin gap. The intuition for this 
is that an elementary magnetic excitation can be viewed as 
"breaking" a valence bond by replacing it with a triplet (or 
a pair of 5 = 1/2 "spinons"), which costs an energy of or- 
der J if the bond is between nearby spins but a vanishingly 
small amount if the bond is very long. An experimental con- 
sequence of having a spin gap A is that the magnetic sus- 
ceptibility should vanish exponentially at low temperatures: 
X~e-*/ T . 

These issues have gained additional prominence in light of 
recent experiments probing the magnetic properties of the re- 
cently discovered compound ZnCu3(OH)gCl2, also known as 
herbertsmithite. In this material, the magnetic properties are 
determined by the 5=1/2 copper ions which arrange them- 
selves in nearly perfect, widely separated, kagome planes 16 . 
Measurements of the powder magnetic susceptibility at high 
temperatures, when fitted to a Curie- Weiss law, reveal an an- 
tiferromagnetic exchange constant J ~ 200 K 17 . A variety 
of different measurement s 17 ' 18 ' 19 ' 20 confirm that the material 
shows no evidence of magnetic order down to temperatures 
as low as 50 mK ~ 10~ 4 J, which suggests a non-magnetic 
ground state. 

The simplest model consistent with these facts is the two- 
dimensional 5 = 1/2 kagome Heisenberg antiferromagnet 
(KHAF). Indeed, while less is known about the KHAF than 
its counterparts on other lattices, such as the triangular^ it is 
widely believed that its ground state is non-magnetic 21 . Exact 
diagonalization studies indicate a small spin gap22 A ~ J/20 
and below this scale, the spectrum shows a large number of 
singlet states, the number growing exponentially with system 
sizo^ 3 .. These facts suggest a picture where the low energy 
physics of the KHAF, and hence ZnCu3(OH)gCl2, is domi- 
nated by short-range valence bonds and valence bond solids 
are, in fact, among the proposed ground states i 24 ' 25 ' 26 

However, the coherence of this perspective is disturbed by 
the puzzling fact that the material shows no sign of a spin gap 
and the powder magnetic susceptibility does not go to zero as 
T — > 0. While the original experiments 17 , in fact, showed 
the susceptibility continuing to increase even at temperatures 
small compared to J/20, subsequent experiments have sug- 
gested that the susceptibility might eventually saturate^ or 
perhaps decrease before eventually saturating at a nonzero 
value i 2 - 7 - 
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A number of proposals have been made to resolve this dis- 
crepancy including suggestions that the KHAF ground state 
may actually be a gapless liquid state involving long-range 
valence bonds 28 ; the gaplessness is a disorder effect involv- 
ing magnetic defects 29 and/or nonmagnetic impurities 29 ^; the 
true spin Hamiltonian of the material is closer to an Ising 
model 31 ; and that Dzyaloshinskii-Moriya (DM) interactions 
play an important rol e 32 ' 33 . In this paper, we focus on this last 
idea but will comment on the other suggestions further below. 

Microscopically, the DM interaction between spins arises 
due to the spin-orbit coupling. Originally proposed on the 
basis of symmetr y 34,35 , it can be derived microscopically as 
a linear (in the spin-orbit coupling) correction to the standard 
superexchange mechanism 36 . The interaction has the form: 

#DM = £ H bu = £ D« ■ (Si x S,-) (1) 

{ij) (ij) 

where the sum is over pairs of spins on a lattice. As with the 
Heisenberg interaction, the dominant contribution comes from 
the pairs of nearest neighbors, to which the sum is commonly 
restricted. Because the interaction is antisymmetric, we also 
need to choose a convention for how the pair (ij) is oriented, 
i.e. whether it appears in the Hamiltonian as S, x Sj or as 
Sj x Si (see Fig.[T). Eq. (Q]i can be viewed more formally as 
the antisymmetric part of the most general bilinear interaction 
between spins ^ 

The vectors {Dij } are constrained by crystal symmetries to 
follow certain rules. One rule is that if spins i and j both lie in 
a mirror plane of the lattice, then Djj must be perpendicular to 
this plane. For a literally two-dimensional crystal embedded 
in three-dimensions, the lattice itself is a mirror plane so only 
the out-of-plane component D z will be present. However, in 
ZnCu3(OH)gCl2, this symmetry is broken by the (OH) groups 
which mediate the superexchange between Cu ions, so the in- 
plane component D p will also be present. Another rule is that 
= if the midpoint of the line connecting spins i and j 
is a center of inversion. For the square and triangular lattices, 
the midpoint of every such line is a center of inversion so for 
perfect lattices, DM will not be present. However, the mid- 
points of the bonds forming a kagome lattice are not centers 
of inversion 37 so a DM interaction between nearest-neighbor 
spins can exist. The relations between different 's are de- 
termined by the requirement that Dy transforms like a vector 
under a symmetry operation. For the kagome lattice, these 
relations are summarized in Fig.Q] 

The suggestion of DM coupling as a means of reconciling 
the presence of a spin gap in the KHAF with its apparent 
absence in the real material was first explored by Rigol and 
Singh in Refs. Il32ll33ll . The authors noted that the DM cou- 
pling would mix the singlet and triplet sectors so total spin 
is no longer a good quantum number. Instead, the ground 
state will contain 5 = 1 components which contribute to a 
nonzero susceptibility at T = 0. The authors studied the 
KHAF augmented by the nearest-neighbor DM interaction 
given by Eq. ((T). Fitting their results to the experiment, they 
estimated the values of |D*| and \D Z \ to be - 0.2 J - 0.3 J 
and ~ 0.1 J respectively. However, they also pointed out that 
the numerical techniques used in their study were most re- 




(a) (b) 



FIG. 1 : The relations between the directions of the (a) out-of-plane 
and (b) in-plane components of the DM vectors for nearest-neighbor 
interactions in the kagome lattice are summarized. The directions are 
determined by the physical requirement that Hdm preserves the sym- 
metries of the lattice, but also by the convention of how we orient the 
links as they appear in the Hamiltonian. The figure uses the conven- 
tion of Ref. 13311 where the DM interaction on the link (ij) appears in 
the Hamiltonian with the lower spin or, on horizontal links, the left 
spin, being the first member of the cross product. A different choice, 
which is sometimes convenient, is where the links are oriented clock- 
wise around each triangle. 37 In this convention, all of the D z 's come 
with the same sign. 



liable at temperatures larger than 0.3 J. On the other hand, 
one may argue that the observed susceptibility behaves truly 
"anomalously" only at temperatures smaller than the putative 
spin gap, since only at these temperatures do we expect \ to 
start decreasing due to singlet formation. For the KHAF, such 
an energy scale is ~ J/20, 22 Therefore, as Rigol and Singh 
noted in their conclusion, a theory taking into account DM 
interactions in the low temperature regime is much needed. 

In this paper, we are attempting to fill this gap by studying 
i?DM at the T = limit, which emphasizes the quantum na- 
ture of the system. Our main conclusion is that the picture of 
ZnCu3(OH)gCl2 having a low temperature phase dominated 
by short-range valence bonds can be reconciled with the sus- 
ceptibility measurements if we include the effect of a pertur- 
batively small DM interaction. 

Our approach is not specific to the kagome lattice and we 
begin in section|ll]by discussing an analytical method for cal- 
culating the susceptibility of a valence bond system in the 
presence of weak DM interactions on an arbitrary lattice pro- 
vided that the unperturbed Hamiltonian satisfies certain as- 
sumptions. The assumptions we require are: (1) the unper- 
turbed Hamiltonian has a narrow band of low energy short- 
range valence bond states separated from magnetic states by a 
spin gap, (2) that i/oM connects this narrow 5 = band to a 
narrow band of 5 = 1 magnetic states, and (3) that the unper- 
turbed Hamiltonian conserves the total spin. For such models, 
we show that the DM interaction leads to a nonzero value of 
the zero temperature susceptibility. 

Then, in section [Till we apply our approach to a model 
of the spin-dimer compound SrCu2(B03)2, which lives on a 
Shastry-Sutherland lattice. The reason for considering this 
model first is that the assumptions and approximations of sec- 
tion [II] hold exactly there. Then, in section [TVl we consider 
a generalized Klein model on the checkerboard lattic e) 10 ' 38 . 



In this case, we have a Hamiltonian whose ground states are 
known to be short-range valence bond states; whose spectrum 
has approximately the form we assume; and for which the 
overlap expansion is believed to converge fairly well. This 
will provide a level of confirmation that our methods can yield 
a correct order of magnitude in a "real" problem. In section 
|VJ we return to the kagome lattice and ZnCu3(OH)gCl2. In 
section I VII we discuss our results in light of recent numeri- 
cal work, other experiments, and alternative theoretical view- 
points. 



□ 
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(b) 



H. MODEL AND FORMALISM 

In this section, we calculate the effect of a small DM inter- 
action on the low temperature powder susceptibility of a lat- 
tice spin system whose low energy properties are determined 
by short-range valence bonds. Our calculation applies for a 
system that satisfies a few rather general assumptions, which 
will be stated at the beginning of section MB I For simplic- 
ity, we concentrate on the case where the valence bonds are 
always between nearest neighbor spins but the argument is 
more general. We begin by discussing some technical aspects 
of working in a (nearest-neighbor) valence bond basis. 

A. Dimer basis and overlap expansions 

A configuration where every spin is in a singlet with one 
of its neighbors can be represented pictorially as a dimer cov- 
ering of the lattice, as shown in Fig. for a square lattice. 
Therefore, we will refer to these basic configurations {|d Q )} 
as "dimer coverings". The wave functions which will interest 
us can be written as superpositions of these dimer coverings: 

|<£) = Ys^K) (2) 

a 

but this requires some care because the dimer coverings are 
not orthogonal. One consequence is that the inner product of 
\(f)) with a state \ip) = ^2 a b a \d a ) is not simply J2 a a a^a but: 

a,/3 

where fl a p = {d a \dp) is the overlap matrix. 

The magnitude of the overlap of two dimer coverings is 
most easily calculated by overlaying the configurations to 
form their transition graph, as shown in Fig. [2J3. As shown 
in the figure, the resulting picture contains double bonds and 
loops of various even lengths. The magnitude of the corre- 
sponding overlap matrix element is: 

I "a/3 I =2 N 'l[x L > (4) 

i 

where Ni is the number of loops in the transition graph; the 
product is over these loops, L, being the length of the ith 
loop; and x = l/y/2. Thus, the overlap between two ar- 
bitrary dimer coverings will often be small, as the transition 



FIG. 2: (Color online) (a) A general dimer covering of the square lat- 
tice. This represents a wave function, which we write as \d a ), where 
a dimer across a link, shown as a thick red (black) line, means that the 
two spins are in a singlet bond, (b) Overlaying two different dimer 
coverings, shown by the red and blue (black and gray) lines, gives 
the transition graph, which will contain both double bonds, where 
the two coverings coincide, and closed loops of various even lengths. 

graph will contain many long loops, though (for a finite sys- 
tem) never zero. However, there is also a notion of a maximal 
overlap which occurs between two dimer coverings that differ 
by only one minimal length loop. 

This latter observation is the basis of the overlap expan- 
sion, which is an approximation scheme based on treating x 
in Eq. (0]l as a small parameter. 5 For example, the overlap ma- 
trix of a set of square lattice dimer coverings, to leading order 
in this expansion, is: 

fi^T « <W - 2x 4 D Q/3 + . . . (5) 

where \3 a p equals unless the dimer coverings \d a ) and \dp) 
differ by exactly one minimal length loop, which on the square 
lattice has length four. In this case, O a p = ±1> where the 
sign depends on the sign convention we take for the singlets 
(i. e. our choice in whether a singlet between spins i and j is 
written as ^(ifii - HJt) or "^OVl - Hit))- lt can be 
shown that for any lattice, we can choose the sign convention 
so that the overlap of two dimer coverings differing by exactly 
one minimal loop always comes with a negative signJS For 
the square lattice, this means that the entries of the matrix □ 
are always or 1 . We will assume this sign convention in this 
paper. 

The overlap expansion can be applied to a general operator 
O. For example, on a square lattice: 

O a0 = (d a \0\dp) m A a 5 a p - 2x 4 B a0 n afi + ... (6) 

where A a — O aa and B a p are constants. Eq. © is especially 
convenient when A a and B a p are independent of a, (3, 

We will use the overlap expansion in the calculations that 
follow, but we remind the reader that since x is actually 1 / v2, 
the approximation is a poorly controlled one. Nonetheless, the 
expansion has proven to be a useful guiding principle in the 
construction of effective models of valence bond dominated 
phases i 5 i 9 i 39 i 40 We expect that for lattices where the minimal 
loop is large, such as the length six loop of the kagome, the 
approximation should work fairly well at leading order. 
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B. Perturbation theory for susceptibility 

We now discuss our assumptions about the unperturbed sys- 
tem (i.e. the Hamiltonian in the absence of DM interactions). 
Then, armed with the technical background of the previous 
section, we shall derive a perturbation theory for the suscepti- 
bility in the presence of a small DM interaction. 




Hdm will mix the singlet and triplet sectors, we are assuming 
that the excited states {|e 7 )}, for which the matrix elements 
(rf Q |i/DM|e 7 ) ^ 0, are specifically S = 1 eigenstates. The 
simplest way to meet this requirement is via our third assump- 
tion: that the unperturbed Hamiltonian conserves total spin, 
which both the Klein and Heisenberg models do. We also as- 
sume these states form a narrow band. One way this condition 
may arise is if the magnetic excitations can be viewed as local 
disturbances of the low energy states. The simplest example 
of such a disturbance would be to break a single dimer by ex- 
citing the singlet to a triplet. 

The precise way in which each of these assumptions are 
used will be made explicit during the derivation. 



(a) 



(b) 



(c) 



2. Derivation 



FIG. 3: The generality of our first assumption allows us to consider 
different possibilities for the spectrum of the unperturbed Hamilto- 
nian. (a) The low energy sector is comprised of a narrow band of 
states which are superpositions of dimer coverings, separated by a 
spin gap A from the excited states, (b) The same spectrum as in 
(a), but allowing a number of low-lying singlet states to exist in the 
gap. (c) Our second assumption: We assume Hdm connects the nar- 
row band of low energy S = states with a narrow band of higher 
energy S = 1 states. 



1. Assumptions 

Our first assumption is that the spectrum of the unperturbed 
Hamiltonian has a low energy structure resembling Fig. [3^. 
The figure depicts a set of low lying states each of 

which can be written as a superposition of dimer coverings: 



(7) 



These states form a narrow band of width 8 and are separated 
from magnetic excitations by a spin gap A ^ i5. While hav- 
ing a spin gap is crucial to the analysis, our formalism can be 
adapted to a situation where higher energy singlet states occur 
in the gapped region, as depicted in Fig. fj}). A well-known 
class of toy models having (or widely believed to have) this 
structure are those of the Klein-AKLT typei 10 i 38 ' 41 ' 42 These 
are models for which dimer coverings, or a subset of them, 
are zero energy ground states and the collection of coverings 
{ | d a ) } forms a spin-gapped degenerate ground state manifold. 
In fact, we expect Fig. [3^b to be a reasonably accurate cari- 
cature of a spin-gapped system even when there is no obvious 
reason for restricting our attention to nearest-neighbor valence 
bonds. We will discuss the suitability of this assumption (and 
the others) for the KHAF in section [V] but we refer the inter- 
ested reader to the summary in section IIF of Ref. Il40ll . 

Our second assumption is that Hum connects our narrow 
band of low energy singlet states {\n)} with a band of S = 1 
excited states {|e 7 )}, whose bandwidth S e is also small com- 
pared to A, as depicted in Fig. fSJ;. While it is well-known that 



Because the DM Hamiltonian (Q]i does not commute with 
the total spin operator, S = £\ S^, the sum being over sites 
in the lattice, deriving an expression for the susceptibility re- 
quires some care. Consider the Hamiltonian: 



H = H o + H DM - g/x B h • S 



(8) 



where Ho is the unperturbed Hamiltonian discussed in the 
previous section, h is the magnetic field, g the g-factor and 
fiB the Bohr magneton. If {E a , \a)} are the eigenvalues and 
eigenstates of H, then the partition function is: 



Z[T,h] = J2 e ~ Ea{h)/kBT 



(9) 



where k B is the Boltzmann constant and T the temperature. 
The magnetization M may be computed in the usual way: 



dF d 



dh, 
z ^ 



k B T dZ 
Z dh a 



dE, 



Oh 



a e -E a /k B T 



(10) 



and, similarly, the susceptibility: 



/tlv= BW L= fl 2 Z 



1 dZ dZ 



dh v \Z dh^dh v Z 2 <9/i M dh v 

1_ ^ ( d 2 E a 
7, 2.^ 



z 
i 



-E a /k B T 



Xdh^dhv 

1 fdE a dE Q 
~7, 2^ 



-E a /k B T 



Z 2 Yl 



dE a 
dh L , 



\ dhfj, dh v 
,-E a /k B T 





dE, 



dh v 



a e -E,3/k B T 



(id 



To compare with experiments, we need the powder suscep- 
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tibility which is given by Xpowder = j(x a 



Xpowde 



1 J2(KE a )e- E ^ T 



3Z 



a 



2 e -E a /k B T 



1 

Z2 



Ene -E a /k B T 



(12) 



where V h = ^ ^ 

If the total spin S commutes with H, the states {|a)} may 
be chosen as simultaneous eigenstates of the two operators 
and the eigenvalues will depend linearly on the magnetic field, 
i.e. E a ~ h-S Q = h-(a|S|a). In this case, the second deriva- 
tive term on the RHS of Eq. (fTTT i will vanish and the remain- 
ing terms will give T x ^ u ~ [{S»S V ) - (5^) (S 1 ')], which is a 
familiar version of the fluctuation-dissipation theorem. How- 
ever, if S does not commute with H, the field dependence 
of the eigenvalues can be more complicated. For example, 
if the leading field dependence is quadratic in h, then in the 
zero field limit, the term in square brackets on the RHS of 
Eq. (fTTI) will vanish and the second derivative term will be all 
that remains*^ 

Our first assumption was that the unperturbed Hamiltonian 
Hq has a spectrum of the form shown in Fig. [3] Because we 
also assumed that Hq commutes with S, the eigenstates of 
Hq can be chosen to describe the system in the presence of 
a magnetic field. In addition, if g^sh is sufficiently small 
compared to A, the low energy spectrum of the Hamiltonian 
Hi = Ho — <?^sh • S will still resemble Fig.[3]and the low- 
est energy eigenstates will still be the collection of S = 
states {\n)} of Eq. <(7j. However, the S ^ bands will split 
into separate bands indexed by the spin component S h along 
the field direction h. In particular, the triplet band {|e 7 )} 

mentioned earlier will split into three bands {\e\ )}, where 
S h = —1,0,1. The spin gaps of the three bands are given by 
A' 5 > = A — g/iBhS h . From Eq. ( flZb . one may verify that 
at low temperatures, the zero field susceptibility will vary as 
Xpowder ~ J'~ 1 e~ A / T which decreases to zero as T — > as 
expected for a spin gapped system. 

Having formally included the magnetic field exactly (at 
least with regard to the low energy spectrum), we now con- 
sider the effect of a small DM interaction, Eq. (|T), on this 
picture. If D = |D| is sufficiently small (a sufficient, though 
not necessary, condition for "small" is if D is small compared 
to the smallest of the three spin gaps A^ s the spectrum of 
the perturbed Hamiltonian H will still resemble Fig.[3]in the 
sense of a set of low lying states separated by a gap, though, 
since H no longer conserves spin, the gap is no longer a "spin" 
gap. To determine the effect on the low temperature suscep- 
tibility, we need to examine how the low energy eigenvalues 
and eigenstates get modified through mixing with 5^0 com- 
ponents. 

To first order in perturbation theory, the states {\n)} be- 



(k\H DM \n) 



Ek 



\k) 



(13) 



where the sum is over all eigenstates of Hi = Hq — g/ih • S 
except |n). The sum will get restricted by our assumption 
that Hum connects the low energy sector {|n}} with a narrow 
band of S — 1 excited states {|e 7 )}, split by the field into 
three separate narrow bands. Moreoever, because the widths 
of these excited bands and of the low energy sector were as- 
sumed to be small compared to the gap, we may approximate 
the denominators by Ek — E n w A^ s \ In this case: 



Ida) 



-EA^E<4 sfc) i^MK>K sh) >) 

S h 7 / 

^2a na \d' a ). 



(14) 



Therefore, in this approximation, the consequence of having 
Hdm present can be visualized in terms of its effect on the 
individual dimer coverings: 



Id' 



Ea^EK 5 ' 1 HHuuKXef )) (15) 

S h 7 



To proceed, we now examine the effect of the operator Hum 
on one of the dimer states. As the operator is a sum over 

pairwise interactions, fl^M ' mere ^ two cases to consider: 
the link (ij) can either be (a) occupied by a dimer or (b) an 
empty bond (see Fig.|4|i. 



J 



FIG. 4: Two different situations for a term in Hdm to operate on: (a) 
a dimer; (b) an empty bond. 

We use the following notation to indicate the spin state of a 
pair of spins on sites i and j: 

(ij) = -/= Or,nj|,n - H,njt,n) (16a) 



Mo = (n,n3l,n + H.njt.n) (16b) 

foil = *T,rJT,n (16c) 

[ij]-! = U.nJi.f. (16d) 
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- the singlet and three triplets associated with the quantization 
axis n where if(j,),„ denotes the S™ = |(— h) eigenstate of 
the operator Sf and the subscripts —1, 0, 1 indicate the value 
of S" + for the pair (ij). The notation for a singlet (ij) 
does not have a superscript because the singlet state is inde- 
pendent of the choice for n. 

The action of Hum on a dimer covering is most easily cal- 
culated with respect to z quantization axis. Using the above 

(ii) 

notation, the action of the term H DM in case (a), when spins i 
and j are in a singlet is: 



while for case (b), when (ij) is an empty bond, the result is: 



D- 
D+ 



{[«]?[#]5 -[«»]?} 



4V2i 
D z 

~ 4i~ 
D 



:{[ei\-r\jk]o-lei\o\jk}-i} 



(18) 



D+ 



{[ek]' {ij) + (ek)[ij] z } 
:{[efc]?(ij) + (efc)[y]?} 
{[ek] z _ 1 (ij) + {ek)[ij]U}- (19) 



4\/2'. 

To avoid clutter, we have dropped the subscript ij on the DM 
coefficients and D ± = D x ± W". Eqs. ([17) and ([18) show 
that the effect of operator on a dimer covering is to pro- 
mote the dimer(s) emanating from sites i and j from singlet(s) 
to triplet(s). As indicated explicitly in Eqs. ( [T7) and ( fT9) , the 

state -ff^M M")' an£ l nence the state floM|^a)> is an eigenstate 
of total spin with 5 = 1. An immediate consequence of this 
is that (n\Hj3M\n), which is a sum of matrix elements of the 
form (d a \HuM\dp), is exactly zero. Therefore, the energy of 
state \n) will modified by Hum only at second order in pertur- 
bation theory. 

The second order correction to the energy of state \n) is 
given by: 



4 2) = E 

K , l 



A (S<*) 
S h a,P,"f 



\(k\H DM \n)\ 2 

E n — Ek 

a na a n (3{d a \H D M\e^ h) ){e^ Sh) \H D M\dp) 

(20) 



The operators •* = ^2 \e^ s ')(e < - s )| are projection op- 
erators that respectively select for the S h = —1, 0, and 1 
components of Hjyn\dp). Therefore, it is useful to rewrite 
Eqs. (TT7T>— (fT9b with respect to the two-spin states of the h 
quantization axis. If: 



then Eq. dTTb becomes: 

[«1S 



and similarly Eqs. ([T8l and ( [T9t become: 



H D $((ei)(jk)) = — {[e{\1\jk]-i [ei] h -i\jk]i} 



D- 



:{[ei]JM-[ei]g[7*]f} 



4a/2' 



:{N-iL?'fc]o-NoM-i} (23) 



--l-{[ e 4 l fe-) + (efc)[y]5} 



4V2i 



{[eA]?(tj) + (eA)[ij]f} 



where 

D h = D h = D*sin0ros<A+D tf Mn0sin<i+D*cos< 



(24) 



sm W cos 0x + sin t/ sin < 



(21) 



(25) 



D ±±: = -D z sin9 + D x (cos cos 4>=fi sin 4>) 

± iD v (cos 4> =F i cos 9 sin</>) (26) 

The _L superscript refers to the direction perpendicular to the 
magnetic field and D ±+ and D^-~ are complex conjugates 
as the notation suggests. From Eqs. d22l and (f24), we see 
that the amplitudes of the S h = component of HDi&\dp) 
are determined by the set of {D^} while the amplitudes of 
the S h = ±1 components are determined by the collections 

The operator HumP^ s 'Hum of Eq. d20b is a sum of 
link terms such as H D ^V^ S ^ H^ . The matrix element 
(d„ | V^ sh =^ H D C $ I dp ) will be proportional to: 

D*D^ = (D o6 .h)(D ar h) (27) 

while (d Q | J ff^ > ^ (5 ' l=±1) ^DM > M/3> will be proportional to: 

Di b -D^ d + + D^+Dj- = 2[D ah • B cd (B ab ■ h) (B cd ■ h)] 

(28) 

The proportionality constant (which might be zero) will de- 
pend on the overlap of states \d a ) and \dp) and the particular 
links (ab) and (cd) involved. 

In order to proceed, we approximate Eq. ( f20b by the leading 
term in its overlap expansion. This means replacing the matrix 

(d a \HuM'P l ' S ^HuM\dfj) by the diagonal term in (the lattice 
appropriate generalization of) Eq. ((S). We expect this approxi- 
mation to be accurate if the mutual overlaps between the dimer 
coverings entering the superposition in Eq. (|7) are small or if 
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the lattice architecture involves large minimal length loops. In 
this approximation, Eq. ( l20t becomes: 



E n ] ~ E A ( S h) E \ a na 



S h 



(ab),(cd) 

+ 0(.t l »"») (29) 



where L m [ n is the length of the minimal loop which can appear 
in a transition graph for the lattice in question and, to the same 
order, ^ Q |a nQ! | 2 w 1 + 0(x Lmin ). We have also used a basic 
property of projection operators: V 2 = V. 

For a dimer covering \d a ), the set of vectors 

{V^ s ^ Him \da)} Uj) can be determined by looking at 
Eqs. (I22l-d24l>. Eq. ( 1291 involves the sum of all possible 
overlaps between pairs of vectors in this set. These pairs can 
be classified into ten distinct types of combinations which 
are listed in Table U with reference to a dimer covering of 
the generic lattice depicted in Fig. [5] The terms that actually 
arise will depend on the connectivity of the lattice under 
consideration. 

The vector P^^Hfiffl \d a ) resembles \d a ) except the 
dimer(s) emanating from sites a and b have been pro- 
moted to triplets. Therefore, the diagonal matrix element 
(d a \H^V^V^H^>\d a ) = except for when links 
(ab) and (cd) involve promoting the same dimer(s) to triplets. 
This corresponds to cases 1, 5, 6, and 8 in Table U The other 
cases in Table U will contribute at higher orders in the overlap 
expansion, which involve off-diagonal matrix elements. 

Using Eqs. (t22l>-(l23ll. ( 1271 and (|28]), Eq.|29] to leading or- 
der in the overlap expansion, becomes: 



a ( 

+ V [^-(D^-h) 2 ]) 

E(2) / 
(D a6 -D ac 
A[a(bc)]e|<i a ) V 

+ 7?[D afc -D ac -(D ab -h)(D ac -h) 



(3) 

0[(ab)(cd)]e\d a 



Dad ' Dbc 



[D ad -D fcc -(D ad -h)(D bc -h)])| 



where r\ 



1 



+ 0(x L '™) (30) 



The first sum in Eq. (1301 . la- 



belled with superscript (0), is over all links in the lattice and 
is a uniform shift that is the same for every state \n) in the 
low energy sector. The three sums, labelled with superscripts 



(1), (2), (3), are taken with reference to a particular dimer 

(2) 

covering \d a ). These terms will therefore contribute to E n 
according to which dimer coverings enter the superposition in 
Eq. (0. The sum labelled (1) is over all links (ab) that contain 
a dimer. The sum labelled (2) is over all triangular plaquettes 
that contain a dimer, where A[a(6c)] means that the dimer 
lives on link (be). The sum labelled (3) is over square pla- 
quettes D[(ab)(cd)] with two dimers on links (ab) and (cd). 
These last two sums will only occur if the lattice in question 
contains triangular and square plaquettes respectively. 

Eq. (l30l gives the field dependence of the energy levels of 
the states in the low energy manifold. We expect these low ly- 
ing states to dominate the thermal averages of Eq. (fTTT i when 
the temperature and magnetic field energy are small compared 
to the gap A due to the Boltzmann factors. If we therefore 
make the approximation of restricting the sums of Eq. QT| 
to the low energy manifold, Eq. [30] tells us that the leading 
field dependence (at low fields) is quadratic in h. Therefore, 
the zero field susceptibility will be determined by the second 
derivative terms of Eq.Qj] which, as mentioned above, vanish 
for systems where total spin is conserved. 

To compare with experiments, we calculate the zero field 
powder susceptibility in the T — ► limit: 



Xpowder 



6A 3 

7El w «Ei a ««i 2 (E 



(i) 



2V (2) 

^— 'A [a 



(6c)]e|d a 

(3) 



(ab)e\d a ) 
D a 6 • D ac 



D, 



[{ab){cd)]£\d c 



, 'Dad ■ D bc 



(31) 



where w n — er En l kBT is the Boltzmann factor of state \n) 
and the primed sum is restricted to the low energy mani- 
fold. The expression can be simplified if we assume the 
magnitude of D is the same on every link and the products 
D ab • D ac = D 2 cos #a and D a d ■ D bc = D 2 cos are the 
same for every triangular and square plaquette respectively. In 
this case, Eq.[3Tlbecomes: 



Xpowder — 



_ N(g^ B ) 2 D 2 fz + 1 



6A 3 



2 cos 6>a 



A 



+2 cos 



N 

(Ng) 
N 



(32) 



where z is the coordination of the lattice, N the number 
of sites, and (N Ai a) = JfE' w « E a \ a na\ 2 N {An)ta where 
iV(A,D), a is the number of triangular (square) plaquettes 
which contain one (two) dimer(s) in dimer covering \d a ). 

Eq. d32l shows explicitly that the low temperature, zero 
field powder susceptibility does not vanish at low tempera- 
tures but instead approaches a constant that depends on D. 
From Eqs. (fTOb and ( f30b , we find that quadratic field depen- 
dence also implies that the low temperature, zero field magne- 
tization will be zero. Therefore, the picture of a short-ranged 
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t 

FIG. 5: Generic 2D arrangement considered for determining the dis- 
tinct combinations of H^, M , as listed in TableU The specific TableU 
cases which arise will depend on the connectivity of the lattice under 
consideration. 
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TABLE I: Calculating the second order (in DM) correction to the 
energy, \ involves computing all possible overlaps of pairs of 
vectors in the set {VH^ \d a )}(ij) for each dimer covering in the 
superposition that defines the unperturbed state n) (Eq. (J7J). There 
are ten distinct cases to consider, indexed in this table by the link op- 
erators involved, in reference to the generic dimer covering in Fig. [5] 
The possibilities are both links involving dimers, (l)-(2); one link in- 
volving a dimer, (3)-(4); and neither link involving a dimer (5)-(10). 
If both links involve a dimer, they can be either (1) the same or (2) 
different dimers. If one link involves a dimer, then the other link is 
an empty bond that either (3) contains or (4) does not contain one of 
the spins of that dimer. If neither link involves a dimer, then the pos- 
sibilities are (5) the two empty links are the same; the two links share 
exactly one spin while the other two spins either (6) form a dimer — 
this is possible if the lattice has triangular plaquettes or (7) do not 
form a dimer; or the two links do not share a spin (8)— (10). In this 
latter case, the two empty links may be connected by (8) two dimers 
— this is possible if the lattice has square plaquettes; (9) exactly one 
dimer; or (10) no dimers. 



valence bond phase with DM interactions is qualitatively con- 
sistent with the combination of a non-vanishing T — > sus- 
ceptibility and lack of magnetic order observed in experiments 
on the herbertsmithite compound. 17 i 18 i 19 i 20 In the following 
sections, we discuss how this approach fares quantitatively 
in the context of various models where our assumptions are 
known or widely believed to be satisfied. 



III. APPLICATION TO THE SHASTRY-SUTHERLAND 
MODEL OF SrCu 2 (B0 3 ) 2 

We now consider a toy spin model defined on the lat- 
tice in Fig. [6^, which was first considered by Shastry and 
Sutherland (SS) in 1981. 43 This model has been revis- 
ited more recentl y 44,45 in light of the spin-dimer compound 
SrCu2 (BC>3)2, where the copper ions form a lattice (FigJSJ?) 




(a) (b) 



FIG. 6: (a) The Shastry-Sutherland (SS) lattice. In its ground state, 
the spins connected by the diagonal bonds form singlets. The arrows 
on the links denote their orientations in terms of the sign convention 
for Hdm, i.e. which spin comes first in the cross product. The ori- 
entation of the diagonal bonds can be taken as going from the lower 
site to the upper site, though this fact will not enter the calculation, 
(b) The structure formed by the copper ions in SrCu2(BOs)2. This 
is topologically equivalent to the SS lattice. Indicated in this figure 
are the directions of the D vectors on the various links, where the 
orientation convention of (a) has been assumed. 

which is topologically equivalent to the one considered by SS. 
From our standpoint, this model satisfies the assumptions of 
section HI1 almost exactly, which is why we discuss it as a first 
application of our formalism. 

A. Model 

In the SS model, the interaction between neighboring spins 
depends on whether the link connecting them is horizontal, 
vertical, or diagonal. The Hamiltonian is: 

Edias ^ — *horiz,vert 
Si-Sj + J'^2 Si-S,-. (33) 

(io) (v) 

where J and J' are positive constants. If J = 0, the model re- 
duces to the square lattice Heisenberg antiferromagnet, which 
is widely believed to have a Neel ordered ground state. If 
J' = 0, the ground state is a product state of singlets on the 
diagonal bonds: |\&) = Y[{if\(^j)- SS showed that this prod- 
uct state, which we may call a valence bond solid, 78 is an exact 
eigenstate of the full Hamiltonian ( 1331 and, in fact, the ground 
state up to a critical value of J'/ J, which they determined to 
be of order unity. When J' = 0, the basic magnetic excita- 
tions involve replacing singlets with triplets and the spin gap 
is exactly equal to J. In Ref. 14511 . it was shown that the exci- 
tations remain nearly localized with a small dispersion when 
J' 0, except for very close to the transition (in fact, it has 
been shown that even singlet excitations are gapped 46 in this 
parameter range). Therefore, in its valence bond phase, the 
SS model has a spectrum like Fig. [3^, where the low energy 
sector now consists of just one state |\&). 

We now augment Eq. ( T33T > with a perturbatively small 
nearest-neighbor DM interaction. For an ideal 2D lattice, we 
only need to consider the component D z . Moreover, for an 
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ideal SS lattice, D = on the diagonal bonds because the 
midpoints of those bonds are centers of inversion symmetry. 46 
However, in the context of SrCu2(BC>3)2, it has been argued 47 
that a slight buckling of the planes removes these symmetries 
and the leading effect is to induce D x and D y terms on the 
diagonal bonds. Therefore, in anticipation of relating to the 
experiment, we consider the DM interactions shown in Fig. [6] 
In adapting Eq. OTb to this model, a few points should be 
noted. First is the fact that Hum couples the ground state to 
two different bands of magnetic excitations (each of which 
gets split into three bands on application of a magnetic field). 
To see this, it is simplest to consider the J' = case first. 
Referring to Table U we see that case 1 terms connect | *f?) 
to the set of magnetic eigenstates where one of the dimers is 
now a triplet. These states have an energy Ai = J above 
the ground state. Cases 5 and 6 connect |\ff) to eigenstates 
where two of the dimers are triplets (case 8 terms, while not 
forbidden by the lattice, do not occur in our considered ground 
state). These states have energy A 2 = 2 J above the ground 
state. Since these two bands are orthogonal to one another, we 
can consider them separately in the calculation. In particular, 
the argument which used the assumption of having a narrow 
band, can be easily generalized to the present case because 
both of these bands are narrow. For instance, one may verify 
that the proper generalization of Eq. ([T5T l is: 




-E-^E^f^M^lef)}. (34) 
s h A 2 7 

where the first sum is over the S = 1, E = Ai = J band 
while the second is over the S = 1, E = A2 = 2 J band. 
These arguments will still apply for the case where J' ^ 
(when we are away from the transition point) though the bands 
will now acquire a small dispersion and the spin gaps will be 
renormalized 4 ^ 6 - (and A 2 will no longer be 2Ai). We will ac- 
count for this in our calculation by using the actual gap values 
Ai, 2 in Eq. d34l 

Similarly, the equation for the powder susceptibility will 
now consist of two separate parts, each of the form in 
Eq. [3T| corresponding to the two bands of excited states. 
For the term associated with the Ai band, D 2 b = 

D 2 ab = N [{D x f + (Dy) 2 ] /4 while the other two 
sums do not contribute. For the term associated with the A 2 
band, £ (0) D 2 ab = 4N d (D z ) 2 = 2N(D Z ) 2 while the sec- 
ond sum will not contribute because in our chosen state, the 
dimers are only on diagonal links. One may verify that each 
triangle contributing to 53(2) Da& ■ D ac contributes the same 
value, D a (, ■ D ac = —(D z ) 2 and the number of these triangles 
is 2Nd = N. Therefore, we arrive at our final result: 

N{g^ B ) 2 f (D*) 2 + (Dy) 2 8(D Z ) 2 \ 

Xpowder- 12 ^ A? + A 3 2 J 

(35) 

A noteworthy point is that because the low energy sector of 
the unperturbed model contains only one state which happens 



to be a dimer covering, higher order terms in the overlap ex- 
pansion will not enter. Therefore, the only approximations in 
Eq. ( |35l l are those inherent in perturbation theory. 



B. Comparison with experiment 

The most striking feature of the temperature dependence 
of the zero field 79 susceptibility of SrCu 2 (BC»3) 2 is a peak, 
which occurs around T = 15 K, followed by a rapid decrease 
to nearly zero as the temperature is further lowered.— The 
most natural interpretation of this rapid decrease is a reduction 
of entropy as the physics becomes increasingly dominated by 
a non-magnetic ground state. ESR measurements 48 revealed 
two sets of triplet excitations with spin gaps Ai = 35 K and 
A 2 = 55 K respectively. 80 

The SS model has had some success as a theoreti- 
cal description of both the ground and excited states of 
SrCu 2 (BC<3) 2 in zero magnetic field. The analysis of Ref. [45] 
determined that the temperature dependence of the suscepti- 
bility was well modelled by Eq. ( 1331 with a value of J'/ J « 
0.68, which is in the valence bond phase but somewhat close 
to the transition point (though still far enough away to have 
appreciable spin and singlet gaps and fairly localized triplet 
excitations 45 *^ 6 .). In addition, it was show n 46 ' 47 ' 49 that by 
including DM interactions, the model could also explain a 
number of features that occurred in the presence of a mag- 
netic field, such as the appearance of uniform and staggered 
magnetization a 44 ' 50 in fields small compared to the spin gap 
scaled 

Eq. d35l ) predicts the effect of DM in the case of zero mag- 
netic field: as T — > 0, the susceptibility will not decay to 
zero but instead approach a constant. To compare with exper- 
iment, we use g sa 2, the measured values for Ai w 35 K and 
A 2 ps 55 K given in Ref. JH], an estimate of D z ss 2 K taken 
from ESR measurements 46 , and the suggestion of Ref. ll47ll 
that D x,y ~ D z . The molar susceptibility is obtained by 
multiplying Eq. ((35) by the factor Na/N, where Na is Avo- 
gadro's number. With these numbers, we obtain the estimate: 
Xpowder(T = 0) « 5 x 1CP 5 emu/mol Cu. The measured sus- 
ceptibility in Fig. 2 of Ref. 14411 does not go to zero and, in 
fact, begins to show a rise at a temperature after reaching a 
minimum value of around 5 x 10~ 4 emu/mol Cu at roughly 
4 K, which the authors attributed to a small (< 1%) concen- 
tration of impurities. 

Therefore, more experiments on substantially cleaner sam- 
ples are needed to test Eq. (l35l >. We reemphasize that the as- 
sumptions of section are expected to hold well for the SS 
model so to the extent that Eq. d33l with DM included is a 
good model of SrCu 2 (BC<3) 2 , we expect our order of magni- 
tude estimate to be reliable. 



IV. APPLICATION TO THE GENERALIZED KLEIN 
MODEL ON THE CHECKERBOARD LATTICE 

We next consider a toy model defined on the checkerboard 
lattice shown in Fig. |7^. The sites of this lattice are the same 
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FIG. 7: (a) The checkerboard lattice with arrows specifying the ori- 
entations of links with DM interactions. On horizontal and vertical 
links, we take D = D z z (i.e. the sign is always positive) while on 
diagonal links D = 0. (b) The minimal loop on this lattice has length 
8. (c) The next minimal loop has length 12. 



as a square lattice but the connectivity is different: on alter- 
nate square plaquettes, opposite corners are also connected as 
nearest-neighbors by diagonal links (note: the crossing point 
of these diagonal links is not an additional site). The lattice 
can be viewed as a 2D projection of the 3D pyrochlore lattice. 

The checkerboard lattice does not currently have a real- 
ization in an actual material. While there are a number 
of materials 51 which form pyrochlore structures, we are not 
aware of a system where the low energy physics is believed to 
be captured by short-range valence bonds. However, there is a 
theoretical motivation for considering this system. The model 
we are about to describe is among the simplest examples of a 
model for which (a) the Hamiltonian and exact ground state(s) 
are known and well-characterized, (b) the assumptions of sec- 
tion Q]] are expected to either hold or be a decent approxima- 
tion, and (c) the overlap expansion is expected to converge 
fairly well, while at the same time the model (d) captures all 
of the complexities of the most general scenario. In this sense, 
we may consider this system as a "laboratory" in which to ex- 
plore the limitations of our formalism. 

The model we consider is a two-dimensional version^ of 
the generalized Klein model, introduced in Ref. [10] (note: a 
model of a similar type was also considered in Ref. 15211 '). The 
Hamiltonian for this model may be written as: 



H 



^2 h p 



(36) 



where the sum is over crisscrossed square plaquettes and: 



given in Refs. [10] and lT38h shows that the ground state mani- 
fold of Eq. ( 1361 ) consists only of dimer coverings of the lattice 
where every crisscrossed plaquette has exactly one dimer and 
superpositions of such coverings. Note that there are dimer 
coverings of the checkerboard lattice that are not included in 
this set. 

In terms of Fig. [3^, the low energy sector of Eq. fl36l l is a de- 
generate manifold. One may consider perturbing Eq. fl36t to 
arrive at a model with a smaller ground state space. It has 
not been proven, but strongly suspected, that the Hamilto- 
nian fl36j > has a spin gap; demonstrations of a spin gap have 
been made, to varying degrees of rigor, in a number of related 
models i 10 ' 42 ' 53 ' 54 Therefore, there are strong reasons to believe 
that Eq. (f36b satisfies the assumptions required by our formal- 
ism. 

Moreover, the loop structure of the ground state space is 
such that there are strong reasons to expect the overlap ex- 
pansion to converge well in a wide number of situations. 82 
Because the minimal loop has length 8 (Fig. [7}}), keeping 
only the diagonal term of the overlap expansion will often 
be a good approximation. Also, the next minimal loop has 
length 12 (Fig. |7j;), instead of 10, so keeping the leading off- 
diagonal term in the expansion will result in a relative error of 
~ 2;(i2-8) ^ yyhjig m man y other lattices, including the 
square, this error will be ~ x 2 . 

The simplest way of introducing a nearest-neighbor DM 
coupling is to consider an ideal 2D lattice, in which case we 
only have D z on the horizontal and vertical links (D = on 
the diagonal links due to inversion symmetry). Because we 
do not have a specific material in mind, we take for simplic- 
ity D z to be the same magnitude and sign on all of the links, 
as per the link orientation convention in Fig. [7^. In adapting 
Eq. QTb to the present case, we note that all four cases con- 
tributing to that general expression will occur here. The first 
sum will give J] (0) D 2 ab = 2N(D Z ) 2 , where 2N is total num- 
ber of non-diagonal links. Similarly, the second sum will give 
D 2 ab = {D z ) 2 (% - AT diag ) where 7V diag is the number of 
diagonal dimers in the coverings under consideration. There 
will be two triangular terms for each diagonal dimer, and for 
each of these terms D a fc-D ac = (D z ) 2 . Similarly, each square 
plaquette having two dimers will contribute —(D z ) 2 . There- 
fore, the analog of Eq. d32b for the present case is: 



h p = (S 2 )(S 2 -2). 



(37) 



where S z 



Si 



S4 is the total spin of the four 



members of plaquette p. h p is an operator that projects the 
wave function onto its component where plaquette p has the 
maximal spin of 2, the two factors in Eq. d37l i annihilating the 
spin and 1 components respectively. Since Eq. (l36b is a sum 
of projection operators, its eigenvalues must be non-negative 
and any state with zero eigenvalue is a ground state. 

A wave function in which crisscrossed plaquette p has a 
dimer on one of its links will be annihilated by h p since the 
total spin of p can then be at most 1 . Therefore, a wave func- 
tion having a dimer on every crisscrossed plaquette will be a 
zero energy ground state of Eq. ( f36b . A counting argument 



Xpowde 



5N(gfi B ) 2 {D 
12A 3 



z\2 



2{N t 



diag/ 



4(iV E 



N 



5N 
0(x 8 ) 



(38) 



where, (• • • ) denotes thermal average as per the discussion 
after Eq. d3lV 

It would be useful to study Eq. (l36*l l numerically to con- 
firm the existence of a spin gap and to examine its low-lying 
5=1 states, in light of the assumptions of our formalism. A 
comparison of an "exact" calculation of Xpowder for (possibly 
a perturbed version of) Eq. d36*i l with the analytical expression 
Eq. d38l ) would provide a useful validation of our approach 
and is natural topic for further study. 
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V. APPLICATION TO THE KAGOME LATTICE AND 
ZnCu 3 (OH) 6 Cl 2 

In the previous sections, we have applied our formal- 
ism to systems with well-characterized Hamiltonians and 
ground states. In the present section, we finally return to 
the original motivation for this work and discuss the effect of 
Dzyaloshinskii-Moriya interactions on the (short-range) va- 
lence bond physics of the kagome lattice in the context of the 
material ZnCu 3 (OH) 6 Cl 2 . 

We begin by assuming the material is well described, to 
leading order, by a spin 1/2 kagome Heisenberg antiferromag- 
net, though the details of our analysis do not depend specifi- 
cally on the KHAF and will apply for any Hamiltonian with 
the assumed spectral properties of Fig. [3] We review some 
known properties of the KHAF which suggest that this model 
satisfies the assumptions of our formalism. We then discuss 
the results of the calculation of Xpowder in light of recent ex- 
periments. We note that while the KHAF assumption under- 
lies most theoretical treatments of this material to date, there 
is not currently a consensus on what model best describes 
herbertsmithite. 31 



A. Model 



of the above mentioned "perfect hexagon" valence bond crys- 
tal state were studied. For that state, it was found that the 
lowest lying triplets have a spin gap of order 0.08 J but a band- 
width of only ~ 0.01J. This is consistent with the assumption 
of a narrow 5 = 1 band required in section Ull 

Therefore, while currently there is no Hamiltonian on the 
kagome lattice for which it can be explicitly shown that the 
assumptions of section Ull are satisfied, the above facts suggest 
that these assumptions are sensible with respect to the KHAF 
and related models. 

In adapting Eq. ( 1321 to the kagome lattice, cases 1 and 5 
(see Table Q) are straightforward (using z=4 for the kagome) 
and case 9 does not occur. For the kagome lattice, each dimer 
is part of a triangular plaquette and each plaquette contributes 
(D z ) 2 + (DP) 2 cos ^ = (D z ) 2 - \{DP) 2 to sum (2) in 
Eq. ( f3TT >. Therefore, 



Xpowde 



N(9li B f 
4A 3 



(2D 2 - (D z ) 2 ) +0(x 6 ). 



(39) 



We can also use our formalism to calculate the anisotropy 
in x, as would be seen in measurements on single crystals: 



X_ 

X" 



2 1 




(40) 



As noted in the Introduction, exact diagonalization stud- 
ies of the KHAF suggest that the ground state is non- 
magnetic^! 2 ^ with a spin gap of order J/20. 22 The gap is filled 
with an exponentially (in the system size) large number of 
non-magnetic states which appear to become a gapless con- 
tinuum in the thermodynamic limit. It was noted in Ref. 112211 
that a gapless singlet continuum and spin gap an order of mag- 
nitude smaller than J could indicate the importance of longer 
than nearest-neighbor valence bonds. However, in Ref. [55], it 
was noted that at least the gapless singlet continuum, could be 
reproduced by considering the KHAF restricted to the nearest- 
neighbor valence bond (dimer) subspace. Indeed, dimer states 
with gapless singlet excitations have been seen in toy models. 5 
In addition, it was noted that the dimer subspace is the sim- 
plest subspace that captures the feature of an exponentially 
growing number of states. 

More recently, the series expansion study of Ref. 12611 ex- 
amined the energies of various dimer coverings of the kagome 
for the KHAF Hamiltonian and found that a particular one, 
namely the "perfect hexagon" state first noted in Ref. Il24ll . 
was optimal though all dimer coverings were very close in 
energy, with a bandwidth (of order J/50) small compared to 
the spin gap. Their estimate for the ground state energy per 
site compared well with exact diagonalization studies of the 
model. Put together, these facts suggest that spectra such as 
Fig. (3^b, where the low energy physics is determined by a 
narrow band of dimer states, are decent caricatures of what 
happens in the KHAF. Further arguments about the suitability 
of restricting attention to the dimer subspace when studying 
low energy properties of the KHAF may be found in section 
IIF of Ref. IHOll. 

In Ref. 115 611 . the elementary triplet (and singlet) excitations 



B. Comparison with experiment 

In comparing our results with experiments on 
ZnCu 3 (OH) 6 Cl 2 , the first point to note is that DM does 
not induce a zero field magnetization, which is consistent 
with the lack of magnetic order according to a number of 
techniques i 17 - 18 - 19 ! 20 

Eq. d39l is an expression for the T — powder suscepti- 
bility. In order to compare with experiment, we assume that 
D z «fl;A~ J/2022 or possible a bit higher - J/10 5 ^ (the 
gap entering the calculation is the one separating the bands 
connected by Hum for which the spin gap is a lower bound); 
and J ~ 170 K^2 

There are two sets of experiments which suggest that the 
susceptibility saturates in the T — » limit. Ofer et al. 18 mea- 
sured Xpowdei- using /iSR in a 2 kG magnetic field in a tem- 
perature range from around 100 mK to around 200 K. They 
observed a monotonic rise in Xpowder as the temperature was 
lowered and the last two data points indicated a saturation 
Value Of Xpowder (T = 0) ~ 15.7 x 10~ 3 emu/mol Cu. Us- 
ing Eq. ( [391 , this value is consistent with D ranging from 
0.03 J to 0.08 J. However, other interpretations of this data 
have involved the Ising anisotropy 31 and/or stressed the role 
of impuritie s 29 ' 30 ' 57 ' 58 . 

The second experiment is a recent 17 NMR study 27 pur- 
ported to measure the intrinsic susceptibility of the kagome 
planes. Because the energy scale of the 6.5 T field used in 
that study is comparable to our lower estimate for A, it is not 
clear that Eq. ( |39l >, the derivation of which assumed that the 
Boltzmann factors of the excited bands were relatively small 
compared to the low energy dimer manifold, will directly ap- 
ply. However, if the actual gap is closer to the higher end of 
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our range ~ J/10, then Eq. ( f39l > might produce a reliable or- 
der of magnitude estimate. The measurements of Ref. rf27tl 
suggest 59 a T — susceptibility per spin of \ = 0.13 which 
implies a molar susceptibility of ~ 1.1 x 10~ 3 emu/mol Cu. 
From Eq. d39l >, this implies that D ranges from 0.008 J to 0.02 
J. 

We caution against taking these estimates too literally be- 
cause in both cases they are based on a relatively small num- 
ber of data points and in both experiments, it is not clear that 
DM interactions are the only factors at play. Also, because 
these two experiments suggest qualitatively different temper- 
ature dependences for %, they can not both be measurements 
of the quantity we are calculating. However, we would like to 
emphasize that these estimates strongly suggest that the mag- 
nitude of the DM interaction may be significantly smaller than 
previously quoted values. We will discuss this point in more 
detail in the next section. 

While our calculation is for the T — susceptibility, we 
can speculate on what happens at a small but nonzero tem- 
perature. Numerical studies of the KHAF^ 2 - suggest that the 
singlet sector has a linear density of states at the very lowest 
energies. An approximate way of accounting for this is for 
a low energy state with energy e above the ground state, the 
gap entering Eq. (f30b will be A — e, if we assume the triplet 
band is still fiat. The sum over n in Eq. ( f3Tb will become an 
integral over a density of states, the range of integration be- 
ing the effective low energy bandwidth 6, In the limit where 
T <C S <C A, the result will be a modification of Eq. d39t : 
X(T) fa x(T = 0)(1 + CT) where C is a positive constant 
that depends on D and A. Therefore, our picture predicts that 
\ should rise with T so in this sense resembles what is seen 
in Ref. 21. However, more experiments and a more refined 
theoretical treatment are clearly warranted. 



VI. DISCUSSION 

The original motivation for this work was the material 
ZnCu3(OH)eCl2, which appears to have a non-magnetic low 
temperature phase, whose nature has been the subject of much 
speculation. At the beginning of this paper, we noted that any 
non-magnetic phase could be viewed as either a short-range 
or long-range valence bond phase. In the previous section, 
we have shown that a short-range valence bond phase with a 
very small DM coupling provides a picture of the low tem- 
perature phase of ZnCu3(OH)eCl2 which reconciles the lack 
of observed magnetic order, the lack of an observed spin gap, 
and known facts about the KHAF. We now discuss this picture 
in the light of more recent experiments and other theories. 

Recently, Zorko et al^S determined the magnitudes of D z 
and D p based on an analysis of high temperature ESR mea- 
surements. They found the best fit was obtained for D z ~ 
0.08 J and |D P | ~ 0.01J. These values were somewhat 
smaller than the estimates of Rigol and Singh 32 , also based 
primarily on high temperature measurements. One possible 
reason for our discrepancy with these estimates is the obser- 
vation of Imai et al. 20 that the OH bonds, which mediate the 
superexchange between the Cu spins, freely rotate about the 



Cu-Cu axis at high temperature but freeze in random orien- 
tations below 50 K. While the effect of this on J would be 
weak, we expect the effect on D to be more significant be- 
cause the direction of D is determined by the position of this 
OH grouppS Therefore, there is strong reason to suspect that 
the strength of the DM interaction experienced by the system 
at low T is different than the value suggested both by the high 
temperature ESR analysis and the numerical fit of the high 
temperature magnetic susceptibility 32 . We emphasize again 
that our estimate of D is based on low T susceptibility data. 
(However, it remains to be seen if an analysis of the low T 
ESR data will lead to a larger or smaller estimate for D.) 

There is another possible reason for the discrepancy with 
Ref. ll60}l . One could speculate that the ESR analysis of 
Ref. B60I1 was based on the assumption that the symmetric 
anisotropic exchange could be neglected in comparison to 
the antisymmetric DM interaction. While this is often a rea- 
sonable approach, given that the former term is quadratic in 
the spin-orbit coupling while the latter is linear, some recent 
magnetization measurements 31 give reason to suspect that the 
symmetric term might also be important in this material. (The 
importance of accounting for such a term has been pointed out 
by Shekhtman et ali&i^ and has received further experimental 
confirmation 63 .) However, in the present context, we again do 
not know whether accounting for this will lead to a larger or 
smaller estimated value of D. 

One reason to believe that our estimates may not be entirely 
off target is the fact that in the presence of sufficiently strong 
D z , the KHAF is expected to order magnetically. For classical 
spins, the presence of such a coupling, no matter how small, 
would favor a long-range three-sublattice 120° magnetic order 
in the x — y plane (with the chirality determined by the sign 
of D z ) 31 . Curiously, a similar result would appear to hold in 
the quantum limit if one were to adopt a proposed algebraic 
spin liquid description of the KHAFi^Mi However, no signs of 
any magnetic ordering in herbertsmithite have been observed 
thus far in contradiction with both the quasi-classical picture, 
even after accounting for the spin- wave corrections 65 , and the 
algebraic spin liquid picture. 

Cepas et al£2- recently studied the KHAF augmented with 
DM via exact diagonalization. They found that the non- 
magnetic phase observed in the D = cas o 21 ' 22 is stable to 
small DM coupling until at D c ~ 0.1J the system under- 
goes a phase transition to a magnetically ordered phase. The 
value of D z obtained in Ref. 116011 is very close to such a mag- 
netic transition while our range of estimates puts the system 
deeper into the non-magnetic phase, making it definitely con- 
sistent with the current lack of any experimentally observed 
magnetic order or phase transition. Relevance to experiments 
aside, our estimate and Ref. M59I1 together form a theoretically 
self-consistent picture because our calculation assumes a non- 
magnetic phase. 

Assuming the picture of a short-range valence bond phase 
is correct, an obvious question would be what is the nature of 
the phase? It is difficult to address this question using our for- 
malism. As mentioned previously, the overlap expansion will 
be more accurate, and hence our calculation more reliable, for 
a valence bond crystal where the fluctuations are weak (such 



13 



as the columnar state on the square lattice) as opposed to a 
liquid or a valence bond solid where fluctuations are strong 
(for example, a plaquette phase 66 ). Phases based on maximiz- 
ing "perfect hexagons", variants of which were suggested by 
Marston and Zeng 24 , Nikolic and Senthil 25 , and more recently 
by Singh and Huse^ 6 -, are candidate states to which our calcu- 
lations might reliably apply. We hope the present work will 
rekindle interest in these and other short-range valence bond 
phases in the context of herbertsmithite. 

A non-magnetic alternative to the above picture is a phase 
based on long-range valence bonds. As these phases are spin- 
gapless by construction, the absence of a spin gap in the ex- 
periment is explained by fiat. One issue with this approach is 
the aforementioned numerical evidence suggesting the KHAF 
has a spin gap. Another issue is that at least one calculation 28 
of the properties of such a phase, a variant of an algebraic 
spin liquid, shows the susceptibility vanishing at low tem- 
peratures as x ~ T, instead of saturating. An already men- 
tioned recent calculation 64 indicates that DM interactions will 
drive this same phase into a magnetically ordered state, which 
has not been seen in experiment so far. A nonzero suscepti- 
bility can, in principle, be obtained within a long-range va- 
lence bond picture by considering phases that break SU(2) 
invariance 67 , by coupling the liquid phase to impuritie s 28,68 , 
or via a state with a spinon Fermi surface 69 . 

To conclude, we would like to reiterate or central result: we 
have shown that Dzyaloshinskii-Moriya interactions can ac- 
count for a nonzero value for the T = magnetic susceptibil- 



ity in a system dominated by short-range valence bonds at low 
temperatures, while the magnetization remains zero. Applied 
to herbertsmithite, we find that such a picture is consistent 
with experiments and our approach may be used to estimate 
the strength of the Dzyaloshinskii-Moriya coupling. There 
are, however, inconsistencies in attempting to reconcile our 
results to the estimates based on the high-temperature data. 
Clearly, further studies, both theoretical and experimental are 
needed to fully clarify these issues. Progress on this front 
would contribute greatly towards understanding this material. 
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The frustration of the classical triangular Heisenberg antiferro- 
magnet is partially relieved by having the spins arrange in a non- 
collinear 120 degree pattern which turns out to closely model the 
ground state of the quantum model 70 , despite initial thoughts that 
it was a spin liquid 2 . In contrast, for the classical kagome Heisen- 
berg model, the ground state is highly degenerate 71 and there is 
no obvious reason we are aware of for why quantum fluctuations 
should favor one of them, though a particular set of them are ap- 
parently favored classically due to "order by disorder'—. 
In fact, a careful microscopic derivation shows that the a full bilin- 
ear interaction between the spins is akin to the familiar Heisenberg 
interaction: H = Y] JSi ■ Sj, where S' has been rotated about 
T>ij by a certain angl o 61 ' 62 ' 73 . While experimentally confirmed^ 3 -, 
and undoubtfully important in discussions of a spontaneous ferro- 
magnetism, this observation is not essential for the lowest order 
calculations of the uniform susceptibility in the absence of mag- 
netic order - the subject of this manuscript. 
Of course, the fluctuation-dissipation theorem will still be true. 
However, the relationship between the uniform susceptibility and 
spin-spin correlation will now involve more than just the equal 
time correlator^ 4 . 

We might view a valence bond solid as "exotic" if it occurs spon- 
taneously due to quantum frustration. However, in the SS model, 
and experimentally in spin-dimer materials 44 -, the valence bond 
solid phase occurs because some links have a stronger interaction 
than others. The discussion in the introduction where we counted 
the valence bond solid as an example of an exotic non-magnetic 
phase, was in reference to the first scenario. 
The experiments we are referring to (Ref. (U) were actually done 
in a magnetometer with a static field of H — 1.0 T, which for 
spin 1/2 corresponds to an energy scale q^bH ~ 1.3 K. The 
temperature range of the experiment is from 1.7 to 400 K, so the 
results can be interpreted by a zero-field theory except for perhaps 
the very lowest temperature points. 

Measurements of the NMR relaxation rate 44 determined the gap 
to the lowest magnetic state to be A i =30 K, consistent with the 
value of Ai = 35 K given in Ref. 1 48] . Fitting the susceptibility 
data gave a slightly lower value of A = 19 K but these and other 
authors 4 ^ have interpreted the ~30 K values as the "spin gap". 
Therefore, we use this value for our estimate. 
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The spin gap scale is ~ A/q^lb ~ 20 T while the onset of a 
uniform magnetization is seen at around 18 T. We point out that 
these interesting field dependencies occur at fields much higher 
than the 1.0 T field of the magnetometer used in Ref. 14411 . the 
results of which we interpret with our zero field calculation. 
One situation where the expansion might not work so well is if the 



wave function were a liquid state of the Rokhsar-Kivelson type. 5 . 
In this case, the error made in keeping only leading terms might be 
quite large because of the cumulative effect of many small terms. 
This point was also made in the discussion section of Ref. 15^1 . 



